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SUMMARY 


Measurements, principally with a hot-wire anemometer, were made in 
fully developed turbulent flow in a 10-inch pipe at speeds of approxi- 
mately 10 and 100 feet per second. Emphasis was placed on turbulence 
and conditions near the wall. The results include relevant mean and 
statistical quantities, such as Reynolds stresses, triple correlations, 
turbulent dissipation, and energy spectra. It is shown that rates of 
turbulent -energy production, dissipation, and diffusion have sharp 
maximums near the edge of the laminar sublayer and that there exist a 
strong movement of kinetic energy away from this point and an equally 
strong movement of pressure energy toward it. Finally it is suggested 
that, from the standpoint of turbulent structure, the field may be 
divided into three regions: (l) Wall proximity where turbulence pro- 

duction, diffusion, and viscous action are all of about equal importance; 
( 2 ) the central region of the pipe where energy diffusion plays the pre- 
dominant role; and ( 3 ) the region between (l) and (2) where the local 
rate of change of turbulent-energy production dominates the energy 
received by diffusive action. 


INTRODUCTION 


The one aspect of tiirbulent shear flow that stands out most promi- 
nently is the transport of stream properties by turbulent motions. The 
transfer process is fundamental, for it not only shapes the mean-flow 
field through momentum transfer but supplies the mechanism by which 
turbulent motions receive energy from the mean-flow field. The well- 
known phenomenological theories were the first attempts to give analyt- 
ical forms for the transfer mechanism by some simple physical consider- 
ations and thereby succeeded in making predictions about the nature of 
the mean-velocity field. Subsequent experiments, however, have clearly 
demonstrated the inadequacies of these theories, and in a systematic 
discussion Batchelor (ref. l) has pointed out the inconsistencies and 
unreal consequences of the assiomptions involved. Recently, Rotta 
(ref. 2) and Tchen (ref- 3) presented more extensive and deeper analyt- 
ical treatments of the nonisotropic turbulence problem, and, althou^ 
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their results show some agreement with the existing experimental 
findings^ they cannot escape some arbitrariness in assuming the nature 
of the transfer mechanism. It is the general opinion among investi- 
gators in this field that more extensive experimental work on the tur- 
bulence mechanism is necessary before a satisfactory analytical formu- 
lation of the problem will be possible. 

In the last few years various experimental investigations were 
carried out in different types of shear flows. It became apparent that 
the dynamic and kinematic processes governing such flows may be quite 
different. The difference between the so-called free turbulent flows 
and flows past solid boundaries has recently been emphasized by the 
phenomenon of intermittency^ first noted in a jet by Corrsin (ref. k) 
and later more clearly recognized and studied in detail in the wake of 
a cylinder by Townsend (ref. 5). It is apparent by now that in flows 
like wakes ^ jets, and those near the free surface of the boundary layers 
the intermittency is present and seems to play a very important role in 
the transfer mechanism, while in pipe and channel flows it is completely 
absent. 

Fairly large amounts of experimental information have been gathered 
about flows with free boundaries (refs. 6, 7, and 8); especially exten- 
sive is the work by Townsend in a turbulent wake (ref. 5)* The recent 
work in two-dimensional channel flow by the present writer (ref. 9) 
brought out some significant features of flows with solid boundaries, 
but the information was far from being complete. One consequence of 
this study was the realization that in order to obtain a complete 
picture, say of the turbulent -energy balance, a knowledge of flow con- 
ditions in close proximity to the wall was of utmost importance. 

Since no similar investigation, with emphasis on the turbulent 
structure, has ever been carried out in fully developed pipe flow, the 
present investigation was undertaken. A pipe also offered the sim- 
plicity of axially symmetric mean flow and a turbulent field nonhomo- 
geneous in one direction only. At the same time it afforded an experi- 
mentally convenient setup for obtaining a shear flow of large scale. 

A diameter of 10 inches was chosen as sufficiently large for the favor- 
able application of hot-wires and for bringing the region near the wall 
within practical reach. To cover the desired range, two working 
Reynolds nimbers were chosen, these being 500,000, based on the diameter 
and an airspeed at the center of 100 feet per second, and 50,000 at a 
speed of 10 feet per second. The higher Reynolds number minimized 
viscous effects while the lower one magnified the extent of the pre- 
dominantly viscous layer near the wall. This turned out to be a fortu- 
nate choice since experimental difficulties and errors encountered in 
one case were usually absent in the other. 

The purpose was to investigate the nature of turbulence and its 
relation to the mean flow, in particular the rates of transfer, 
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diffusion^ and dissipation of energy. To this end measurements were 
made of the relevant mean and statistical quantities, including their 
spatial distributions from within the laminar sublayer to the center 
and energy spectra at several radii. Thus the purpose was to reveal 
details of flow structure that could not be found from studies of mean 
flow alone, and yet stop short of a direct attempt to find the nature 
of the transfer mechanism. It is hoped ^ however, that advances toward 
the latter goal will have been made throu^ the present work. 

The present investigation was conducted under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics- The author wishes to express his gratitude to Dr. G. B. 
Schubauer for his constant interest and encouragement throughout this 
work. The valuable discussions with Dr, C. M. Tchen and Mr. P. S. 
Klebanoff are also much appreciated. Thanks are due to Miss Z. W. Diehl 
and Mr. K. D. Tidstrom who carried out the nujnerical computations and 
preparation of the report. 
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pipe radius, 4.86 in. 

dimensionless turbulent -kinetic -energy diffusion 


rates; Dq - 


aid 


IT 3 r dr 
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D_ = — ; rv 
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voltage fluctuations across hot-wires 

fraction of turbulent energy u^ associated with kj^, 
cm^/sec^ 

fraction of turbulent energy v^ associated with kj^, 
cm^/sec^ 

fraction of turbulent energy w^ associated with k^^, 
cm^/sec^ 
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dimensionless gradient diffusion rate of turbulent 


kinetic energy , 


1 d 


d u^ + 


Ut 


4- r dr dr 


one -dimensional wave number in direction of mean 
flow, cm“^ 

mean pressure at any point in pipe 
mean pressiare at exit of pipe 

dimensionless turbulent -pres sure -energy diffusion 
rates, (PD)„ s i i (PD)^ s ^ 1 ± 


dimensionless turbulent -energy production rates} 

/’t, \ uv dU /_ \ uv dU 

(Pr )o = a — — , (Pr ) = v — r- — 

instantaneous value of pressure fluctmtions 
dynamic pressure at pipe center 

Reynolds number based on diameter of pipe and velocity 
at pipe center 

correlation coefficient of u-f luctuations at two points 
displaced in a radial direction 

coordinate in radial direction; r = 0 corresponds to 
pipe center 


friction distance parameter, r'U^/v 
time 

mean velocity at any point in pipe 
maximum value of mean velocity 
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Ut 

Ui, ^5 

U, V, w 

u' , v' , w' 

V, ¥ 

¥ 0 , ¥t 


p 

friction velocity; Uj- = -v 

total (mean-plus -fluctuating) velocities in x, r, and 
qp directions, respectively 

instantaneous values of velocity fluctuations in 
X, r, and cp directions, respectively 

root-mean-square values of velocity fluctuations in 
X, r, and cp directions, respectively 

mean velocities in radial and azimuthal direction, 
respectively 



dimensionless turbulent -energy dissipation rates; 



where Uj^ refers to three velocity fluctuation 
components u, v, and w and x^s, to three coordi- 
nates; repeated indices indicate siommation 


dimensionless direct-viscous -dissipation rate^ 



X dissipation length parameter 

V kinematic viscosity of air 

7t total (mean -plus -fluctuating) pressure 

p air density 

cp azimuthal coordinate 



IH 
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ANALYTICAL COWSIDERATIOWS 
Equations of Motion of Pipe Flow 


The continuity equation and Reynolds equation in cylindrical 
coordinates for inconqjressible mean flow have the following form: 


^ + i + 1^=0 

r dr r dcp 


( 1 ) 


and 


„ dU „ du ¥ du 1 dp /d 2 id — id — \ 

U — + V — + = I — u + ruv H uw 1 + V V^U 

r dr r dcp 


dx dr r dcp P dx \dx 


,, dV „ dv ¥ dV ¥^ 1 dP /d — 

U — + V— + — rr = UV + 

dx dr r dcp r p \dx r dr 

1 ™ v('v^ - ^ 

dcp r / j,2 J.2 acp / 


id 2 

rv^ + 


> ( 2 ) 


„dW „d¥ ¥d¥ V¥ idP /d — , 

U — + V — + = — UW + 

dx dr r dcp r dcp \dx 


2 dN ¥ 


d — • 1 d 2 2vw\ , i^j , 

— vw + w I + vjVw + 

dr rdcp r/ \ r2dcpr2, 


where 
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In a fully developed turbulent pipe flow the conditions are 
(a) V = 0 and W = 0 

ocp 

(c) The velocity field is independent of the coordinate x 
Equations (2) therefore become 

1 Id — /d% 1 du\ 

= ruv + vl + 

P cbc r dr \dr^ r dr / 


1 ^ A. rv^ + 

P 3r ^ dr r 


0 = - 


d 

dr 


w 


2w 

r 


Integrating the last equation and using the houndary condition w = 0 
at r = a it follows that w = 0 for all values of r- Thus the 
Reynolds equations for the turbulent pipe flow reduce to 


1 ^ 
P 5x 


1 d_ 
r dr 


r^uv 


dr / 


(5a) 


d w£ 

p ^ r dr r 


(5b) 


Differentiating equation (5b) with respect to 

^^/Sr 5x = 0. Thus is independent of r 

and (5b) readily integrate to 


X one finds 
and equations 


(5a) 


2 p 3x 


-r^uv - V + A(x) 


(4a) 


p 


(4b) 


a 


r 


dr + B(x) 
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The boundary conditions are at r = 0 


and at r = a 


uv = 0 


dr 


= 0 


uv = 0 


dr 


v^ = 0 


Also let P = 0 at x = 0 and r = a. Then from equation (4a), 

A(x) = 0, 

1 _ 2 2 

P dx " ‘a 

and, integrating, 

- = Ut^x + C(r) ( 5 ) 

p a 

From equations (4b) and ( 5 ), B(x) = -- U_^x. The equations finally 

a 

■become 


uv 



( 6 a) 



( 6 b) 


It is interesting to note, as already pointed out by Kampe de 
Feriet in the case of flows between parallel walls (ref. 10), that the 
shearing stress uv and the mean velocity occur together in one equa- 
tion and the normal stresses v^ and w^ and the mean pressure occur 
together in the other. From the experimental point of view this has the 
advantage that the equations furnish a method of checking the absolute 
accuracy of the measurements of the Reynolds stresses if the mean- 
velocity and pressure distributions are known. 
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Energy Equations 

The momentim equations may be written in the form 


^ 1 ^ ^ 1^2 1 ^ 1^3 ^ 1^2 ^ _1 ^ 

r Pdx ^ 


bx. br T b(f) r P^X^r^ 



^2,^5 Su2U5 2 ^5^ 2U2U5 2 bji , vfv^, - 2 . i- 

bx. br r bq> / 


Multiplying the three equations by U2, U2, and u^, respectively, one 

gets 


^ 1 ^ ^ 1 ^ 2 ^2^3 ^ ^ 1 ^ - ^ u — 

9 k 9 r r pl>^' 
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2v 


1^2 

2^“" -\cb. 




r 2 \ 9 cp / 


9hi2ii2^ 2 ^^^3 " ^5^) 2 ^ 

+ _ — + + = __ Up — + 

9k 9r r dcp r 


p" 2 ar 


2v 


iv®v ^ 


2 -\ 3 x 



' 9 u 2 \^ 2 /^ 2 \^ 

kSt / r 2 \ 9 cp / r^ 


2 ^3 
r2 ^"2 Sep 


9u2U-2^ ^2^3^ 1 ^3^ 2u 2U^^ 2 ^3 3n 

^ + + iL. 

9k 5r r9cp r pr9cp 


2v 


12 2 
- V u,"^ 
2 5 






9 u 5 


l.{^\ _ ^ + -L u ^ 

/ r^ r^ ^ ^cp 
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Introducing the velocity and pressure perturbations 

U-j^ = U + u 


Xlg = "v 


It = P + p 

averaging, and adding the three eqimtions there is obtained the energy 
equation for the t\nrbulent motion 


— au Id /u^+v^ + w^,p\ 

uv — + rv + ^ = 

dr r dr V 2 P / 


V 

r dr dr 


d u^ + v^ + 


vlh ^ 

r\r 



(8) 


where in the last term Uj^ refers to the three velocity fluctuation 
components and xj, to the three coordinates; the repeated indices 
indicate suinmation. 

The second term on the ri^t-hand side of the equation may be 
neglected compared with the other terms for all values of r, since 

(a) near the center it approaches zero as r > 0 as may be shown 

from the continuity equation; and (b) everywhere else its order of 

magnitude vu^/d^, where D is a characteristic length of the order of 
the pipe radius, is smaller than, for instance, that of the last term 

of the eqtiation; that is, vu^/a^ » vu^/d^, since « 1. 

The equation may now be rewritten as 


uv 



1 

r 


d /u^+v^+w^ p 

_ 4 - _ 

dr \ 2 P 


V 

r dr 


r 


d u^ + v^ + w^ 
dr 




= 0 


2 


+ V 


(9) 
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The energy equation thus obtained has essentially the same form as 
that given by Von Karman for the case of a turbulent flow between two 
parallel plates (ref. 11 ). The first term corresponds to the rate of 
production of tiirbulent energy by the action of the shearing stresses. 
The second term represents the rate of energy change due to transfer of 
both kinetic and potential energies by the radial velocity fluctuations 
and is usually referred to as the diffusion term. The third expression 
may be regarded as a gradient type of energy diffusion and is important 
only very near the wall. The foiorth term expresses the rate of energy 
dissipation into heat by action of the viscosity. 


EQUIPMENT 
Wind Tunnel 


The investigation was carried out in the experimental setup shown 
in figure 1. A centrifugal blower having a capacity of 12,000 cubic 
feet per minute and powered by a I 5 -horsepower constant-speed motor was 
placed inside a large pressure box. The box was octagonal in cross sec- 
tion and was 10 feet hi^ and 16 feet long. The air was discharged from 
the blower into the space surrounding the blower and then passed throu^ 
a honeycomb, three screens, and a large contraction cone. The honeycomb 
consisted of hexagonal cells 2 Inches across and 8 inches long. The 
screens had 2h meshes per inch and a wire diameter of 0.0075 inch. The 
contraction cone, which was made out of hardwood and was circular in 
cross section, reduced the air passage from A 8 Inches to I 8 Inches in 
diameter. The cross section was then further reduced from I 8 Inches to 
10 inches in diameter through a 90° elbow. After the elbow a short 
elastic coupling was used to prevent the transmission of vibration from 
the pressure box to the pipe. In order to minimize any flow irregulari- 
ties due to the elbow, another screen and an 8 -inch-long honeycomb were 
installed in the entrance sections of a 25 -foot-long steel pipe. Fur- 
thermore, in order to accelerate the bo\mdary-layer growth, the pipe 
wall was artificially roughened by gluing floor-sanding paper to the 

surface along a length of 2^ feet. With this arrangement it was found 

from the measured mean-velocity distribution at the end of the steel 
pipe, that is, after an entrance length of about 50 diameters, that the 
flow was fully developed turbulent. Following the steel pipe, a 16-foot- 
long seamless brass tube was attached having an inside diameter of 
9.72 inches. This served as the actual test section. 

The speed of the tunnel was regulated by throttling the intake of 
the blower with adjustable vanes and by venting the pressure box. In 
the early stages of the investigation it was realized that for the very 
low velocity pipe experiments, where the intake of the blower had to be 
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nearly closed^ the heat dissipated in the pressure hox was sufficient 
to affect the flow in the test section. This difficulty was overcome 
hy supplying the air to the system hy a ventilating fan placed on the 
side of the hox as shown in figure 1. 


Hot-Wire Equipment 

The basic hot-wire equipment used d-urlng the Investigation is 
described in detail in reference 12. The hot-wires were made of plat- 
inum drawn by the Wollaston process. Later during the investigation 
platinum-rhodium wires (90 percent platinum and 10 percent rhodium) 
were used and were found to be very satisfactory. The wire diameters 
were generally 0.0001 inch; only when the noise problem was very criti- 
cal were f;tner wires (0. 00005 inch) used. For the measurements of the 
longitudinal components of the velocity fluctuations the length of the 
wires ranged from 0.01 to 0.025 inch. 

Special care was taken in building the X-type of wires for the 
ineasiirements of the cross components v and w. In order to minimize 
wire-length effects and to be able to work very near the wall, the size 
of the wire holder had to be as small as possible. With the use of 
prongs made of fine Jeweler's broaches the dimensions of the holder 
head were cut down to approximately 0.015 inch by 0.005 inch, the wires 
having a length about 0.025 inch. 


Traversing Mechanism 

The traversing mechanism simply consisted of a micrometer screw on 
which the hot-wire support was fastened. The support could be rotated 
in a plane perpendicular to the air flow so that the hot-wire could be 
adjusted parallel to the wall. 

The zero reading of the traversing mechanism (r' = O) was found by 
placing the hot-wire close to the wall (approximately 0.01 inch away) 
and measinring the distance between the wire and its image in the polished 
wall by an ocular micrometer. Since the curvature of the wall was small, 
the space between the wire and the wall could be taken as half the 
observed distance. 


PROCEDURE AND RESULTS 


While fully developed turbulent flow was readily obtained through- 
out the l6-foot test length, tedious and time-consuming adjustments had 
to be made to remove secondary effects and obtain axial symmetry. The 
task was particularly difficult at the low Reynolds numbers where the 
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effect of temperature gradient was most felt. This as well as disturb- 
ances from the elbow had to be eliminated. After the final adjustments, 
pressure and velocity surveys made along the l6-foot test length showed 
that the velocity field at each section was the same and symmetry about 
the axis had been obtained. All the final measurements were made inside 
the tube 2 to 4 inches from the exit. 


Measvirement of Mean Velocity and Pressure 

The experiment was conducted at two Reynolds numbers, 50^000 and 
500 , 000 , corresponding to maximum mean velocities of approximately 10 
and 100 feet per second, these velocities being varied sli^tly in order 
to compensate for the daily changes in air viscosity and density. In 
order to cover the wide range of static and dynamic pressures, an 
inclined manometer 5 feet in length was used. With benzol and a slope 
of 10 to 1 the manometer was sensitive to 0.005 centimeter of water by 
direct reading. With a traveling microscope and separate scale the 
sensitivity was increased to 0.0001 centimeter of water. 

The pressure distribution in the direction of the flow was meas- 
ured throu^ pressure taps located every 2 feet along the brass tube. 

The resixlts are given in figure 2. 

Small total -head and static tubes were made of 0.04— inch -diameter 
nickel tube stock with 0.005 -inch wall thickness. The tip of the total- 
head tube was flattened to an opening of 0.006 inch. The static tube 
was placed approximately O.ii- inch above the total-head tube. This 
arrangement was used for measuring the conqjlete velocity distribution 
in the low Reynolds nimber case. In the vicinity of the wall a correc- 
tion had to be made because of the large level of velocity fluctuations. 
The correction, having the form 


“corr = - 


+ V^ + W^ 


U. 1 - 


U*- 




was of the order of 5 percent or less. 

For the hi^ Reynolds number a hot-wire was used to explore the 
mean-velocity distribution near the wall. Here again a correction for 
velocity fluctuations had to be made because of the nonlinear behavior 
of the hot-wire. This was accomplished by an approximate graphical 
method using the known static response curve of the wire (voltage 
against velocity) and the known root-mean-square value of the voltage 
fluctuations. The maximum correction was about 10 percent, the correct 
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mean velocity being hi^er than the observed. At the low Reynolds nm- 
ber the velocity profile determined by the hot-wire was unreasonably 
low in the vicinity of the wall after applying the correction. Except 
for this^ the agreement between the hot-wire results and those obtained 
with the total-head tube was good. 

The measurements are presented in figures 5 and 1<-. In figure If 
the dashed lines indicate the wall velocity gradients computed from the 
pressure drop. The agreement with the measured values is satisfactory. 


Measurements of Turbulence Levels and Shearing Stress 

The three components of the velocity fluctuations u', v', and w' 

and the turbulent shearing stress uv were obtained by standard tech- 
niques described in reference 6. 

Figures 5 to 8 give the detailed results. In figure 8 the solid 
lines represent distributions calculated by equation (6a) using the 
independently measured mean -velocity gradients and u' and v’. It is 
seen that for the high Reynolds number the agreement between these and 
the directly measfired points is very good while for the low Reynolds 
nfjmber the measured points are somewhat higher. It should be mentioned 
that lower acciiracy of all the measurements in the low Reynolds number 
flow is to be expected mainly because of the difficulty of forming con- 
sistent time averages because of the inherent low-frequency fluctuations. 


Measurements of Triple and Qxiadruple 
Velocity Correlations 

Measurements of the triple and quadruple velocity correlations are 
given in figures 9? 11. The basic technique described in refer- 

ence 12 was essentially adopted for these measfxrements . The schematic 
diagram of the electronic circuit is given below: 
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The two hot-wire signals ej^ and 62 sre fed into two identical 
compensated amplifiers after which e^^ is squared. The preliminary 

squaring circuit uses the nonlinear characteristics of two balanced tri- 
odes as employed in Townsend's circuit (ref. 15 ). The output of the 
mixing circuit, which is a simple resistance network (ref. 12 ), gives 
simultaneously the sum and difference of the inputs. Squaring cir- 
cuits 1 and 2 contain a series of diodes properly biased to give the 
square of the inputs. If the signals before the mixing circuit are 
equalized the ratio meter reads directly the correlation coefficient 

It is seen that in order to obtain the triple velocity correlations 
by this method the mean foiirth power is necessary. This could have been 
avoided, of course, if instead of using the ratio meter the outputs of 
the two squaring circuits had been recorded separately and then sub- 
tracted. It was felt, however, that the method adopted gave more con- 
sistent results. The mean fourth power was obtained simply by feeding 
the amplified hot-wire signal into a squaring circuit and then squaring 
the output again by reading the final output on a thermocouple meter. 

Unfortunately, for the lower Reynolds number some inconsistency in 
the measurements was found using the above technique. In comparing the 
values of the double -correlation coefficient uv/u'v' obtained by this 
setup (without the preliminary squaring circuit) with those obtained by 
the conventional method, they were foimd to be from 20 percent to 50 per- 
cent lower everywhere except near the wall. At the high Reynolds nimiber 
the agreement was very satisfactory. Although, because of time limita- 
tions, it was not possible to trace definitely the cause of this dis- 
crepancy, it is believed to be due to a difference in the low-frequency 
phase shift of the two amplifiers. Consequently, the large difference 
between the triple -correlation distributions at the two Reynolds num- 
bers in figure 11 should be considered to be due to experimental error. 


Measurements of Various Dissipation Terms 


The expression for the rate of energy dissipation as it appears in 
the turbulent -energy equation (9) has the form 



I 1 ^ 1 

r2 \Bcp ) \Scp/ r^ \Scp / 


( 10 ) 
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The first three terms were obtained by the differentiation method intro- 
duced hy Townsend (ref. 15 ). By making the assumption 


A = A 

~ " U 3t 


(11) 


he writes 


/^\2 /:v \2 

I A. 1 - -L.1 A 1 

\Sbc/ u2\^/ 


Thus hy electronically differentiating the hot-wire signal e oc u (or 

V, or w), (■^j easily obtained. A recent paper by Lin 

(ref. 14)' shows that this assumption is valid if no mean-velocity gra- 

p 

dient exists and (u/U) « 1. In a shear flow he gives an additional 

condition 


U 


3x 


Su 


In the present measurements these conditions were foimd to be satisfied 
with the exception of a region inside the laminar sublayer. Also during 
the course of a boundary-layer investigation at the UBS an experimental 
verification of the validity of this method was made at O.O5 times the 
boimdary-layer thickness. 


The fourth and seventh terms of equation (10 ) were 
tially by a method first suggested by Taylor (ref. I5), 


obtained es sen- 
in which 


„ _ uu(r) ^ 

Rj. _ 1 

u'u' (r ) 



as r — > 0, where u and u(r) denote fluctuations at distance r 
apart. Thus by measiiring the correlation coefficient Rj, for small 

/N \2 

values of r, \^/ calculated. During the present investi- 

gation the accuracy of the method was greatly improved by adopting the 
following technique: The above equation may be rewritten as 


2 

u2 ^ u^(r) ~ ^(^)] 2. _ ^ 

2u'u'(r) 2u'u'(r) 2u^'^'^ 
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Since ; for small values of v, 



one has 



Several values of a can he easily measured hy placing two hot-wires 

at various small distances r apart. Then le calculated 

_2 

from the slope of a strai^t line in a plot of a against r . 


The remaining four terms of equation (lO) could also have been 
obtained with the method described above. However, the wire arrange- 
ment necessary for such a technique made its application impractical. 
Comparing the first three terms it was found that they satisfy the 
isotropic relations fairly well except near the wall (figs. 12 and I 3 ). 
The other two measured terms, while not too different from the first 


three in the center region, are considerably higher 


^especially 



near the wall. This, of coiorse, is not surprising since dissipation 
lengths or microscales in the radial direction are expected to be 
smaller because of the presence of the wall. For want of a better 
procedure it was assumed that mean-square derivatives with respect to 
a given direction separately satisfied the isotropic relations 


2 




/ V 2 / \2 / v2 

1 _ 2 f^\ _ 1 

r^ r^ r^V^T/ 


It is not possible to estimate the accuracy of this assumption at 
particular points. However, one can determine the error made in the 
total energy dissipation of the pipe cross section. Integrating the 
energy equation (9) across the pipe there is obtained 


2 rt 

o2. 


ruv 


dr 


dr = 


2it 


a 


Wr dr 


( 12 ) 
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which states that the total energy produced at a given section is dis- 
sipated at the same section because of the homogeneity of the field in 
the x-direction. Equation (12) was found to be satisfied within 10 per- 
cent. It is remarked that if the isotropic relation between all terms 
had been assumed and the dissipation calculated in the usual way using 


w = I5v(|y 

the right-hand side of equation (12) would have been smaller hy approxi- 
mately a factor of 2.5* It is still an open question^ however, how 
accurate the approximation is from point to point. 


Energy-Spectrum Measurements 

The amplified hot-wire signal was fed into a Hewlett-Packard wave 
analyzer with a frequency range of 10 to l6,000 cycles per second. The 
analyzer selectivity characteristics were obtained by calibration; two 
fixed band widths were chosen, their effective values being 11 and 
k2 cycles per second approximately. The output of the analyzer was fed 
into a thermocouple circuit. In order to make the measurements inde- 
pendent of the amplifier frequency response the thermocouple output 
readings of the hot-wire signal were matched by a known sine-wave 
input . 

The measurements are presented in figures 14, 15^ and l6. Wire- 
length corrections were applied only to the u* -spectra, using the method 
described in reference l6. 


DISCUSSION 


As already pointed out in the ” Introduction, " the present under- 
standing of the turbulence phenomenon is not complete enough to be able 
to attack directly the problem concerning the natiu’e of the tirrbulence 
transfer mechanism that establishes a stable, nondecaying turbulence 
field such as the pipe flow.. The principal aim of the present work was 
an attempt to obtain an over -all pictxire of the turbulence structure 
without tiying to understand the detailed mechanism responsible for this 
structure. To this end the energy equations of both the mean and tur- 
bulent flow serve as a useful guide. The former in effect expresses 
the relative magnitude of the mean-flow energy loss to the turbulent 
field as compared with losses due to dissipation by direct molecular 
action. The latter gives a relation between the different forms of 
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turbulent -energy rates such as rates of production, dissipation, and 
diffusion. Unfortunately it is not possible to obtain from these equa- 
tions any explicit information on how these various energy changes take 
place. However, it is hoped that once their relative importance in 
different regions of the turbulent field is better understood one mi^t 
attack the ultimate problem of the turbulence mechanism with greater 
confidence . 

A different, but already well-established, method of approach was 
also tried. Instead of investigating the different rates of changes of 
the txrrbulent energy at a given point of the field, its spectral distri- 
bution in the wave-number space was examined. The merits and short- 
comings of this method will be discussed later. 


Mean-Energy Balance 

Multiplying the integrated momentum equation ( 6 a) by the mean- 
velocity gradient dU/dr, the energy equation of the mean flow is 
obtained: 



2 ^ 
dr 




Thus the energy available because of the pressure drop along the 
pipe is partly converted into turbulent energy and is partly directly 
dissipated. The two terms on the rl^t-hand side of the equation are 
compared in a nondimens ional form near the wall in figure I7. There 
are two points that should be noted here: (a) The bulk of the direct 

viscous dissipation takes place in a very narrow region, r* < 15; 

(b) the position where the laminar shearing stress is equal to the tur- 
bulent shearing stress (viscous dissipation equal to turbulence pro- 
duction) is found to be approximately at the same point where the maxi- 
mum amount of energy is produced (r* = II.5). This point is usually 
referred to as the edge of the laminar sublayer. It is seen that not 
only is the bulk of the energy taken from the mean flow directly dissi- 
pated but a considerable portion of the total turbulence production 
(20 percent) also takes place here. 

It is quite apparent from this picture that, in order to obtain "the 
complete picture of the turbulent -energy balance, conditions near and 
within the laminar sublayer have to be known. This, of course, is a 
very difficult task from the experimental point of view. 
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Turbulent-Energy Balance 

From the measurements presented in the previous section all the 
terms of the turbulent -energy equation (eq. (9)) can be calculated with 
the exception of the pres sure -energy diffusion. Unfortunately, since 
this term is small everywhere except near the wall, its determination 
from equation (9) is very inacc\irate. Nevertheless, it is possible to 
obtain in a qualitative way an over-all picture of the complete 
turbulent -energy balance at a cross section. Since the experimental 
problems and errors are different in the two Reynolds number flows they 
will be discussed separately. 

Low Reynolds nixmber flow . - In order to be able to study the flow 
conditions very close to the wall it is necessary to produce as thick 
a viscous layer as possible. This may be done by carrying out the meas- 
urements at very low mean speeds. Certain compromises, however, have to 
be made. First it is much more difficult to establish good flow condi- 
tions at low speeds, and second it is more difficult to carry out time- 
averaging processes in the measurements of statistical quantities. On 
the other hand the absolute magnitude of the dissipation terms may be 
better established since problems such as wire-length effects, amplifier 
frequency response, and noise are not so critical. 

Figiu’es 18 and I9 give the distribution of all the energy terms in 
the low Reynolds number flow. It should be noted that in forming dimen- 
sionless quantities for the coordinates, the characteristic length a 
was chosen in figure I8 and v/Ut was used in figure I9 representing 
conditions near the wall. The correspondence between respective coordi- 
nates is indicated parenthetically in figure 19- 

The following interesting conclusions may be drawn: 

(1) Throu^out the whole cross section, with the exception of the 
center region, the rate of energy production at a point is approxi- 
mately balanced by the rate of energy dissipation. 

(2) All the various energy rates reach a sharp maximum near the 
edge of the laminar sublayer. 

(3 ) This edge appears to be also the region from which the txirbu- 
lent kinetic energy is diffusing both toward the pipe center and toward 
the wall and toward which the pressure energy is transported. 

There is some question about the direction of the pressure diffu- 
sion in the center region of the pipe. It was mentioned in the previous 
section that the measured triple velocity correlations are believed to 
be too low in this region. Assuming that their dimensionless value is 
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approximately the same as the values found in the high Reynolds nTuriber 
flow - an assumption which holds approximately true for all other meas- 
ured dimensionless statistical quantities - kinetic-energy diffusion is 
estimated as shown with a dashed line in figure l8. It is believed that 
the corresponding pressure-diffusion distribution (dashed line) comes 
closer to the actual picture than that obtained from the directly meas- 
ured points (solid line). This would Indicate that the direction of the 
pressure diffusion is toward the laminar sublayer and its value near 
the center is close to zero. 

Hl^ Reynolds number flow . - As already pointed out, accurate meas- 
urements of the various dissipation terms in the case of hi^ Reynolds 
number flow become much more difficult, since, because of the extent of 
the energy spectrum to high frequencies, amplifier -tube noise and wire- 
length effects become increasingly more critical. Although all precau- 
tions were ftiade to minimize these effects, the calculated rate of energy 
dissipation is believed to be too small. From the energy-spectrum meas- 
urements it was possible to infer that considerable dissipation takes 
place at frequencies as hi^ as 50>000 cycles per second where the 
response of the compensated amplifier ceases to be linear. From the 
information gathered in the low Reynolds number flow, an attempt was 
made to estimate the dissipation. The following assumptions were made: 
(l) In the vicinity of the wall where similarity with respect to 
Reynolds number was foimd (this will' be discussed in detail later) the 
dissipation can be obtained directly from the low Reynolds nvimber meas- 
urements; (2) the errors in percentage are the same across the pipe. 

This was found to be fairly closely true when the dissipation measure- 
ments were repeated using various hi^-frequency cut-off filters. With 
these assimiptions and with equation (12) the estimated dissipation was 
obtained. Figures 20 and 21 show the directly measured and estimated 
values. The figures also indicate that the picture of the energy bal- 
ance is similar to that obtained in the low Reynolds number flow. 


Energy-Spectrum Considerations 

In the. study of isotropic turbulence the concept of the energy spec- 
trum proved to be a useful and a convenient one both from an analytical 
and a physical point of view. By considering the distribution of the 
turbulent energy in the wave-number space it is possible to give in many 
instances a simple pictirre of some basic physical processes taking place 
in the field. An example is the commonly accepted picture of the energy 
transfer from smaller to larger wave numbers. 

The utility of the spectrum function is much more restricted in a 
nonisotropic turbulence field where it cannot be contracted into a scalar, 
as in the isotropic case, but has to be treated as a tensor quantity. A 
restrictive factor in all cases, and particularly in nonlsotropic tur- 
bulence, is the fact that measurement is possible of only the so-called 
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Taylor, or time, spectrim. This means that, even if Taylor's hypothesis 
(eq. ( 11 )) is fulfilled, only a spatial energy distribution over a 
siirface = Constant can be measured. This, of coiarse, prohibits 

detailed knowledge of the point -by-point energy distribution in the 
wave -number space. Nevertheless, it is possible to draw some conclu- 
sions from this type of measurements. 

The measured u' -spectra (fig. l4) indicate that the spectral distri- 
butions obtained not too close to the wall manifest similar behavior 
over a wide range of wave numbers. Specifically, they vary as the 
- 5/5 power of the wave niamber k^^ over a considerable range. This is 

clearly indicated in figure 22 where the function is shown 

to be a constant over the range 1 < ki^ < 2k. The two sets of points 

shown in this figure indicate the magnitude of the length corrections 
applied to the directly measinred values. This, of course, is the same 
type of behavior as that of the equilibriimi range of the energy spectrum 
in an isotropic field first predicted by Kolmogoroff (ref. I 7 ). One may 
infer therefore that for sufficiently large turbulent Reynolds nijmber 
flows (u'>/v > 200 for all measurements) and for the mean-velocity gra- 
dients not too large there exists a wide range of wave numbers in which 
the energy represented by the u' -component is transferred from smaller 
to larger wave niJimbers without being significantly influenced by the 
turbulent -energy production mechanism or by viscous dissipation. This 
apparently is true, even though local isotropy does not exist in this 
range as will be seen later. It is further seen that this range varies 
with r' /a and at r' /a = 0 . 0082 , where the mean- velocity gradient is 
already large, the spectriun shows a different distribution. It might 
be mentioned that in this case there is a rather wide wave-number range 

where the spectrum varies closely as kp”^ as predicted by Tchen 
(ref. 3 )- 


The conclusions one may draw from the measurements of the v'- and 
w' -spectra are significantly different. In order to clarify this dif- 
ference a comparison is made between the v' -spectrum measiired at the 
center of the pipe and that calculated from the measured u'-spectr\jm 
using isotropic relations (fig. 25). It is seen that there is a large 
energy deficiency in the low wave -number range, while the energy content 
of the higher wave numbers is much larger. Furthermore, the -5/5-power- 
law type of distribution is completely missing. Similar statements may 
be made for the v' -spectra (and w' -spectra) at points other than the 
center of the pipe. An explanation of such behavior is, of course, dif- 
ficult mainly because the one -dimensional nature of the measurements 
conceals the detailed spatial energy distribution. However, it is hard 
to conceive that the geometry of the field alone would be responsible 
for such a behavior. Figure 23 may suggest that wave numbers receive a 
large portion of their v' energy not only by means of the usual trans- 
fer mechanism from smaller to larger wave mambers. By considering 
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the energy equations for the turbulent components separately^ the term 
most likely to influence the v' -spectrum in the wave-number range con- 
sidered is the pressure term transferring energy from one velocity com- 
ponent to the other. Its form p suggests that it is associated 

with wave numbers larger than those where the bulk of the tixrbulent 
energy is produced and smaller than the ones where viscous dissipation 
is important. This term' therefore may very well be responsible for the 
larger v' energy present in the wave-number range in question. It is 
difficult to draw any conclusion about the role played by the diffusion 
terms; in general they are thou^t of as representing a low-frequency 
phenomenon and therefore would Influence the v' -spectrum in the low 
wave -number range only. 

It should be emphasized that the above discussion is merely specu- 
lative in nature and a considerable amount of experimental work espe- 
cially in connection with the presstire terms is necessary to be able to 
give a quantitative picture of the energy balance spectrumwise. 


General Considerations 

On the basis of the presented measurements the flow field in the 
pipe may be divided into three regions exhibiting different behavior 
from the point of view of turbulence structure. 

V7all-proximlty range . - Jfeasurements near the wall indicate that the 
well-known wall -proximity law of Prandtl (ref. l8) for the mean velocity 
may be extended to the fluctuating-velocity field also. Thus by using 
and v/U- 5 - as the characteristic velocity and length parameters, the 

various velocity distributions become independent of the Reynolds n\mi- 
ber in the approximate range 0 < r* < 50 (figs. 2h to 26 ). 

The various energy rates as they appear in the equation of tm’- 
bulent energy are found to be of equal relative importance; their magni- 
tudes are much larger than those in other regions of the turbulent field 
and they therefore play a dominant role in the energy balance over the 
entire field. 

Center portion of pipe . - Since the direct effect of viscosity on 
the turbulent field is negligible in the center portion of the pipe, it 
may be expected that the distributions of the fluctuations expressed 
relative to the characteristic velocity are independent of the 

Reynolds number. Figures J, and 8 show this to be the case. 

This region is further characterized by the fact that it receives 
a large portion of its turbulent energy by diffusive action. 
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Intermediate region .- The intermediate range extends from r* 100 
to r’/a « 0.1. By considering the turbulent -energy balance in this 
region (fig. 21 ) it is seen that the rate of energy diffusion is much 
smaller than that of production or dissipation. This means that the 
energy produced here is locally dissipated. It should be noted that 
this is one of the implied assumptions of the mixing- length theories. 
Another assumption, namely the existence of statistical similarity 
between the velocity coB 5 )onents, is, however, difficult to accept in 
view of the discussion in the section Energy-Spectrum Considerations. 
Thus, although the gradient type of momentum transfer involved in mixing- 
length theories has more experimental support in this region than in 
other portions of the field, its use is not con 5 )letely justifiable. 


SUMMARY OF RESULTS 


Fully developed flow in a large pipe was found to provide a very 
useful medium in which to study the structure of turbulence in shear ■ 
flow. Embodied in the term "structure" are the interactions between 
turbulent motions and mean flow, and the various transfers and movements 
of energy from point to point and from mean flow throu^ the spectrian of 
turbulent motions. 

The following are the major results: 

1. The importance of a detailed knowledge of conditions in the 
close proximity of the wall was demonstrated. 

2. Using the similarity parameters Ut and v/Uj- the flow field in 
this region was shown to be independent of the Reynolds number. 

5 . The various turbulent -energy rates, such as production, dif- 
fusion, and dissipation, were found to reach a sharp maximum at the 
edge of the laminar sublayer (r* 12), their magnitude being of equal 

relative inq)ortance but much larger than those in other regions of the 
pipe cross section. 

4 . It was found that there exist a strong transfer of kinetic 
energy away from the edge of the laminar sublayer and an equally strong 
movement of pressure energy toward it. 

5 . In the center region of the pipe, the characteristic length and 
velocity parameters were shown to be a and Ut. 

6. In the region of large mean-velocity gradients but outside of 
the dissipative region (between r* 100 and r' /a O.l) energy dif- 
fusion was found to be small compared with the turbulent-energy production. 
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7. The spectimm measurements indicated that in regions where the 
mean- velocity gradients are not too large the u' -spectra vary as the 
- 5/3 power of the wave number over a considerable wave-number range. 

8. A considerably different distribution of the v'- and w' -spectra 
was measured and an explanation of such a behavior was attempted. 


National Bureau of Standards, 

Washington, D. C., October 28, 1952. 
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Figure 1.- Schematic diagram of test setup. 
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Figure 3.- Mean-velocity distribution. 




Figure k.- Mean- velocity distribution near wall. Dashed lines are computed 

from pres sure “drop measurements. 
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Figure 7.- v' and w' distributions 
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Figure 11.- Triple-velocity-correlation distributions. 
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Figure l6.- w' -spectra. R = 500^000. 
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Figure l8.- Turbulent-energy balance. R = 50^000. Dashed lines are 

estimated. 




Figure 19." Turbulent-energy balance near wall. R = 50^000* Numbers in 
parentheses along ordinate correspond to subscript o, as in figure 18. 
Dashed lines are estimated. 
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Figure 20.- Turbiilent -energy balance. R = 500,000. Dashed lines are 

estimated. 


Figure' 21 .- Turbxilent-energy balance near wall. R = 500^000* Dashed lines 

are estimated. 
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Figure 25*- Comparison of measured and calculated v'-spectrtim. R = 500,000 

r’/a- = 1-0. 
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Figure 24-.- Mean-velocity distribution near wall. 
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Figure 25* - u*, and w* distributions near wall. 
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Figure 26.- Tixrbulent shearing stress and double-correlation coefficient 

near wall . 
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